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Measurements of the phonon density of states by inelastic neutron time-of- flight scattering and 
specific heat measurements along with first principles calculations, provide compelling evidence 
for the existence of an Einstein oscillator {"rattler") at ujei ~ 5.0 meV in the filled skutterudite 
Ybo.2Co4Sbi2. Multiple dispersionless modes in the measured density of states of Ybo.2Co4Sbi2 at 
intermediate transfer energies (14 meV < cj < 20 meV) are exhibited in both the experimental and 
theoretical density-of -states of the Yb-fiUed specimen. A peak at 12.4 meV is shown to coincide 
with a second Einstein mode at uje2 ~ 12.8 meV obtained from heat capacity data. The local modes 
at intermediate transfer energies are attributed to altered properties of the host CoSba cage as a 
result of Yb-filling. It is suggested that these modes are owed to a complementary mechanism for 
the scattering of heat-carrying phonons in addition to the mode observed at uiei ~ 5.0 meV. Our 
observations offer a plausible explanation for the significantly-higher dimensionless figures of merit 
of filled skutterudites, compared to their parent compounds. 

PACS numbers: 63.20.dd, 63.20.Pw, 61.05. fg, 84.60.Rd 
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I. INTRODUCTION 

Skutterudites, RM4Sbi2 (R = rare earth element; M 
= Fe or Co), are an important class of materials for stud- 
ies of fundamental physics problems, since some of them 
exhibit unconventional superconductivity,"'^ while others 
show heavy fermion behaviors.^ In the applied area, cer- 
tain types of filled skutterudites, such as iron and cobalt 
antimonides, hold the promise of becoming excellent can- 
didates for real-life thermoelectric applications, particu- 
larly in the 400 - 800 K range. In the mid-90's Slack pro- 
posed the phonon glass-electron crystal (PGEC) model, '^ 
which established the lattice thermal conductivity, k^, 
as an important tuning parameter in the quest for ther- 
moelectric materials with high dimensionless figures of 
merit, ZT.^ The filled skutterudites have been shown to 
have much lower kl's than their parent compounds, since 
heavy atoms, such as La, Ce. Yb or Tl fill the cages of 
the open host structure and are expected to exhibit lo- 
calized atomic vibrations, called "rattlers", which scatter 
heat-carrying acoustic phonons. The presence of such lo- 
calized oscillations superimposed on collective motions of 
plane wave modes are among a few known examples of 
Einstein modes in nature. The "rattler mode" scenario 
seemed to emerge as the correct model for explaining 
experimental neutron time-of-flight scattering data and 
specific heat measurements in the filled skutterudites,^"® 
thus providing a plausible explanation for the much lower 
lattice thermal conductivities in the latter compared to 
their parent compounds. 

However, the findings of a recent inelastic neutron 



experiment on Ce and La-filled skutterudites suggested 
that the guest eigenmodes are coupled to lattice degrees 
of freedom.^ At the same time, some very recent theo- 
retical calculations by Veithen and Ghosez^° have sug- 
gested a possible coupling of guest and host dynamics 
at low transfer energies in the densities of states in La- 
filled skutterudites, consistent with the findings of Koza 
et al..^ Furthermore, the same authors have suggested 
that heavier elements, such as Tl, may possess single Ein- 
stein modes in the density of states, while arguing that 
what has been interpreted as a second guest oscillator 
(i.e. Qe2 «150-200 K),^^*^ could be due to coupled lat- 
tice vibrations. ^"^ Therefore, understanding the physical 
basis of the second Einstein mode, Qe2, is relevant for 
understanding the large differences in the ZT's between 
filled skutterudites and their parent compounds. 

In this article, we aim to contribute to the understand- 
ing of filled skutterudite thermoelectrics by comparing 
the experimentally-measured density of states (DOS) of 
Ybo.2Co4Sbi2 from time-of-flight inelastic neutron scat- 
tering data with a theoretical Density Functional Per- 
turbation Theory (DFPT) calculation of the DOS in 
Ybo.25Co4Sbi2, and correlating the Einstein tempera- 
tures obtained from a heat capacity measurement on 
Ybo.2Co4Sbi2 with local modes observed in both the the- 
oretical and experimental DOS's. 

The paper is organized as follows. In Sec. II we dis- 
cuss the materials preparation methodology, review the 
neutron time-of-flight instrumentation and data analysis 
utilized in the derivation of the neutron-weighted den- 
sity of states of Ybo.2Co4Sbi2, the theoretical approach 



for the derivation of the density-functional calculation of 
the DOS in Ybo.25Co4Sbi2, and finally, the measurement 
technique of the heat capacity of Ybo.2Co4Sbi2 with an 
outline of the procedure for obtaining the Einstein tem- 
peratures of the studied sample. The main findings of 
the joint experimental and theoretical work are given in 
Sec. III. Finally, we conclude with a summary in Sec. 
IV. 



II. MATERIALS AND METHODS 
A. Materials Synthesis 

A polycrystalline Ybo.2Co4Sbi2 sample was prepared 
via a solid state method. High-purity Yb (pieces 99.9%, 
Alfa Aesar), Co (Cobalt Slut 99.95 %, Alfa Aesar) and 
Sb (shots, 99.9999%, Alfa Aesar) were mixed in stoichio- 
metric quantities. The mixture was sealed in a carbon- 
coated quartz tube under 300 mbar Ar pressure. The 
quartz ampoule was heated to 1000°C at 2.5°C/min, then 
slowly heated to 1100° (at w l°C/min) and left for about 
30 h. Subsequently, the ampoule with the homogeneous 
molten liquid was removed from the furnace at 1100°C 
and quenched into a water bath. The ingot which formed 
as a result of the quench was placed in a furnace and an- 
nealed at 680°C for 168 h in order to form the correct 
crystallographic phase. Subsequently, it was powdered 
and spark-plasma sintered into bulk pellets at 50 MPa 
and 620° C for 2 minutes in a protective Ar atmosphere. 
The high density (> 99.5% theoretical density) pellets 
consist of single phase Ybo.2Co4Sbi2 compound, as con- 
firmed by X-ray diffraction and transmission electron mi- 
croscopy analysis. 



B. Inelastic Neutron Time-of-Flight Measurements 
and Analysis 

The inelastic neutron time- of- flight neutron scattering 
studies were performed on the Cold Neutron Chopper 
Spectrometer (CNCS) of the Spallation Neutron Source 
(SNS) at Oak Ridge National Laboratory. The sample 
(w 9.15 g of polycystalline millimeter-sized pieces) was 
loaded into an aluminum housing fastened to a suspended 
aluminum frame. A cadmium mask was placed in front 
of the holder exposing only the sample to the incident 
beam in order to reduce the background signal. Prior 
to the scattering event, the neutrons were selected by 
a high-speed double disk chopper with sensitive phase 
adjustment operated in the high-resolution mode of the 
instrument. Subsequently, the scattered neutrons were 
collected by detectors placed 3.5 m from the sample. 
Neutrons with incident energies of 9.09 meV provided 
for a usable range from -5 to approximately 45 meV on 
the neutron energy gain side. The energy resolution at 
zero energy transfer was 0.186 meV at FWHM. The sam- 
ple and the empty aluminum can were each measured at 



ambient temperature (w 300 K) for 8 h in order to ob- 
tain a satisfactory signal-to-noise ratio in the measured 
intensities. 

The experimental results were analyzed as follows: The 
recorded neutron counts were normalized with respect to 
proton charge on target and sorted in 1 fiscc bins. Sub- 
sequently, the dynamic structure function S{q,uj) was 

o -1 

integrated in the g-range from 3 — 6 A . The back- 
ground spectrum, comprising of scattering through an 
empty aluminum holder, S^{uj), was subtracted from the 
sample spectrum, S{uj)^ and converted into a precursor 
to the density of states as follows: 



GP(a;) 



[S (w) - S'' (w)] • w 
coth (^) - 1 



(1) 



, where 6*^(0;) and w stand for precursor to the density of 
states (in arbitrary units) and the neutron energy gain (in 
meV), respectively, while /3 = (fcBT)"\" GP{uj) is the 
neutron-weighted density of states and it deviates from 
the true density of states G{uj) for a number of reasons. 
Most notably, multiphonon scattering contributions are 
also embedded in G'p(w), which need to be removed in 
order to obtain the true density of states G{oj)}^^^^ 

A number of techniques for removing multiphonon 
scattering from Gp{uj) have been employed. ^'^ These 
range from the simplest removal of a constant, 
to treating the multiphonon part in the incoherent 
approximation. ^■^^^^ Here, we calculated the multiphonon 
scattering to all orders in the incoherent approximation 
using a phonon expansion procedure described in de- 
tail elsewhere. ^^ Not knowing the phonon DOS a pri- 
ori^ the phonon expansion was calculated using a the- 
oretical partial phonon DOS calculated for each atom 
type from first principles. The multiphonon scattering 
contributes about 2% relative to the 1-phonon scatter- 
ing and was mostly flat. Since the multiphonon contri- 
bution was small and since the final phonon DOS was 
fairly similar to the theoretical phonon DOS, no further 
iterations were needed. In order to remove the residual 
background signal from air scattering, a third order poly- 
nomial was fitted through the low-energy (0 — 4 meV) 
and high energy (37.55 — 44.95 meV) tails of G'p{uj) and 
subsequently subtracted from G^(a;) in order to obtain 
G[uj). The density of states is expected to vary as w^ at 
the lowest transfer energies < cj < 3 meV, and vanish 
above the high-energy optic modes (io > 36 meV), pro- 
viding the boundary conditions for the fit. Finally, G{uj) 
was normalized and the results displayed in Fig. 1(a). 

The experimental density of states (see Fig. 1(a)) was 
used for the derivation of the sample specific heat by 
performing the discrete sum over the energy exchange, 
LOi, according to: 
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FIG. 1. (color online). Shown is the remarkable correspon- 
dence between (a) the G(iu) of Ybo.2Co4Sbi2 derived from 
the inelastic neutron time-of-flight neutron scattering experi- 
ment and (b) the results of the density functional perturbation 
calculation of the phonon DOS in Ybo.25Co4Sbi2 (plotted in 
gray). Individual contributions of the different elements are 
as shown: Sb - red, Yb - green and Co - blue (b). The ex- 
perimental data in (a) exhibit modes which account for those 
predicted by the theory. A clearly-defined mode is visible at 
(a) u ~ 5 meV, agreeing with (b) a distinct peak at uj « 5.52 
meV in the theoretical G(cj). 



, where h and Aw stand for Planck's constant and the 
discrete energy separation between data points (Aw = 
0.1 meV), respectively. The heat capacity derived from 
the DOS {G (uj)) offers a plausible consistency check on 
the results of the direct measurement. The results of the 
conversion are exhibited in Fig. 2. 



C. Theoretical Approach 

The "Frozen Phonon" method^ ^ was adopted in or- 
der to obtain the lattice dynamic properties of Yb-fillcd 
skutterudites. We started from a totally-optimized cell of 
Ybo.25Co4Sbi2 (2x2x2 supercell of Co4Sbi2 together 
with 2 Yb atoms) by using rigorous computational pa- 
rameters: lxlO~^ cV for self-consistent calculation and 

o 

1x10^ cV/A for Hcllmann-Feynman force convergence. 
Then, the atoms were given small finite displacements 

o 

C/q = ±0.01 A, and the number of displacements was de- 
termined by the symmetry of the system. All the energy 
and force calculations were carried out under the projec- 
tor augmented wave (PAW) method. -'^^'^^ implemented in 
the Vienna ab initio simulation package (VASP).^° 

The force constant matrix elements are the second 
derivatives of the total energies with respect to the pair- 
wise finite displacements. The dynamic matrix D(q), 
where q is a phonon wave vector in the first Brillouin 
zone, was derived by taking the Fourier transform of the 
force constant matrix. The square of the phonon fre- 
quency w(q,z^) and the eigenvector e(q,t/) were computed 
by solving the following equation: 



D(q) • e(q, v) = uj^{q, v) e(q, v) 



(3) 



, where v is the index of vibrational modes. Then, the 
partial phonon density of states was calculated by: 






(4) 



, where n, i and /i indicate the number of q points, atom 
number and fraction of thrce-dircctional parts, respec- 
tively. In our calculation for Ybo.25Co4Sbi2, we used a 
Gauss broadening function g(uj) with a smearing factor 
of 0.1. The resulting DOS is shown in Fig. 1(b). 



D. in-situ Heat Capacity Studies of Ybo.2Co4Sbi2 

The heat capacity, Cp, of Ybo.2Co4Sbi2 was measured 
on a Physical Properties Measurement System unit by 
Quantum Design, utilizing its Heat Capacity Option. A 
small piece of material was removed from the same ex- 
perimental sample for the inelastic neutron time-of-flight 
measurement, which was subsequently cut and polished 
to a parallelepiped {m = 0.0244 g; 2.195 mmx 1.917 
mmxO.931 mm). 

In order to assess the relative contributions of lattice 
phonons and "rattler" Yb atoms, we measured the heat 
capacity of the Yb-filled sample with respect to temper- 
ature. The focus was predominantly on the low temper- 
ature range of the data (6 K < T < 53.5 K), such that 
the experiment would probe low-energy heat-carrying 
phonons without the obscuring influences of high-energy 
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FIG. 2. (color online). The heat capacity, Cp, of 
Ybo.2Co4Sbi2 divided by the temperature, T, plotted as a 
function of 7* (i) measured m-situ (orange •) and [ii) de- 
rived utilizing Eq. 2 from the neutron scattering DOS (red 
+). Inset: the low temperature section of the heat capacity 
of the sample and the relevant temperature range for a fit 
with one Debye and two Einstein modes (from the direct in- 
situ measurement). The black line is the fit through the data 
points. 



fix) =j + P-x + A- Qei^ ■ x-(3/2) • 



„©Bi/\/3; 



Q&E2/V^ 
'e©E2/V^ _ I) 



(5) 



, where x represents 2^, 7 is the electronic contribution 
to the specific heat, /3 = ^^^^^''b . (e^)-3 ^jth N - 
the number of host atoms comprising the skutterudite 
system, and &d — the Debye temperature. &ei and 
Qe2 stand for the Einstein temperatures of the two guest 
modes necessary for the fit with spectral weights A and 
B, respectively. The results of the fit are examplified by 
the black line through the data points in the inset of 
Figure 2. 

The quality of the fit is significantly better with the 
inclusion of the second Einstein oscillator: x^ = 1-2 x 
10-3 with only e^i, and x^ = 2.7 x 10"^ with both 
Qei and Qe2 (see Table I). The contributions of all those 
parameters to Cp/T are graphically displayed in Figure 
3. 

The same series of fits with the above-mentioned mod- 
els have been performed on the heat capacity derived 
form the experimental DOS (red '-I-' plot in Fig. 2). 



phonons and anharmonic effects. ^^ The heat capacity, 
Cp, of Ybo.2Co4Sbi2 was divided by the temperature, 
T, and plotted as a function of T^. The results are dis- 
played in Figure 2 in J/K^-mole of atoms (plotted with 
orange circular signs). 



A number of different approaches have been employed 
in the analysis of heat capacity data of filled skutterudites 
in the quest for a "rattler" eingenmode.^"^ In a number 
of them, the Cp of the parent compound is subtracted 
from the one of the filled skutterudite and the resultant 
curve is fitted with the expected contributions of one or 
two Einstein modes. ^'* These approaches to analyzing 
heat capacity data in filled skutterudites are based on 
the assumption that the only difference between parent 
and filled skutterudites is the one attributed to a rattling 
mode, but fail to examine other possibilities, such as hy- 
bridizations between the vibrational motions of the guest 
and antimony atoms forming the cage.^°'^^. We followed 
the analytic approach used by Mandrus et al.^ and fitted 
the Cp/T of the experimental sample Ybo.2Co4Sbi2 in 
the range 6 K < T < 53.5 K directly with a function of 
the form:^ 
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FIG. 3. (color online). The relative contributions of the elec- 
tronic specific heat 7 (in pink), the Debye term, /3 (in purple), 
and the two Einstein modes, 0_ei (in red) and Qe2 (in green), 
obtained from the fit through Cp/T versus T^ of the direct 
in-sjiu-obtained data with Eq. 4. The total Cp/T is shown in 
cyan diamonds. All the data are plotted with their respective 
errors. If not visible, the error is smaller than the diameter 
of the data points. 



Type of Fit 


One Debye Mode Only 


One-Debye-One-Einstein 


One-Debye-Two-Einstein 


7 (J/mol-K^) 


0.0133 ± 0.0015 


0.0054 ± 0.0003 


0.0042 ± 5.2x10"'' 


/3 (J/mol-K") 


6.85x10"'' ± 1.7x10"'' 


3.15x10"'' ± T.TxlO"' 


1.52x10"'' ± 3.6x10"'' 


A (J/mol-K) 


- 


5.356 ± 0.134 


1.622 ± 0.041 


Obi (K) 


- 


103.66 ± 0.89 


70.04 ± 0.47 


B (J/mol-K) 


- 


- 


8.914 ± 0.104 


6^2 (K) 


- 


- 


149.02 ± 0.90 


X' 


0.049286 


0.0011986 


2.6706x10"'' 


R 


0.94963 


0.99881 


0.99997 



TABLE I. Shown are the results of three different fits on the heat capacity measurement of Ybo.2C04Sbi2 from the in-situ data 
(orange • and line): (1) with only one Debye mode, (2) with one Debye and one Einstein modes, and (3) with a Debye and 
two Einstein modes, as shown in Fig. 4. R stands for the correlation coefficient of the fit. The Debye temperature, Qd was 
derived from /? (Od — 307.65 ± 2.39 K), and it is comparable to the one obtained in Lao.gFe3CoSbi2.^ 



III. RESULTS AND DISCUSSION 

The in-situ studies of the specific heat in Ybo.2Co4Sbi2 
(see Fig. 2 orange • data) focused on the low- 
temperature domain (6 K < T < 53.5 K) in order to 
avoid the obscuring influences of anharmonic cfFccts and 
interference from high-energy optic modes. ^^ The data 
was fit in the specified regime with (i) just one Debye 
mode, (ii) one Debye and one Einstein modes and (iii) 
one Debye and two Einstein modes, in order to show the 
goodnesses-of-fit for all these cases (see Table I). 

Similarly, the heat capacity data obtained from the 
conversion of the experimental DOS (Eq. 2) was also 
fitted with Eq. 5 (one Debye and two Einstein modes) 
in the 6 K < T" < 54 K range, in order to check the 
consistency of the data in Table I. The following results 
were obtained: 7' = 0.0169 ± 0.0004 J/mol-K^; p = 
1.44 X 10-5 ± 5.7 x lO"'^ J/mol-K''; A = 1.440 ± 0.055 
J/mol-K; e^i = 60.616 ± 1.020 K; B = 6.309 ± 0.156 
J/mol-K; Qe2 = 146.53 ± 2.020 K; x^ = 7.3886 x 10"^ 
and R = 0.99995. By comparison, the fit of the same data 
with just one Einstein mode resulted in x^ = 0.00035394 
and R = 0.99776. 

The in-situ heat capacity measurement is suscepti- 
ble to error introduced by parasitic heat losses, black- 
body radiation effects, while the DOS conversion is sen- 
sitive to background subtraction and the resolution of the 
CNCS. However, the remarkable agreement between the 
fitting coefficients obtained from the fits with two Ein- 
stein modes on both the in-situ and DOS-derived heat ca- 
pacity data, compared to those of any other model used, 
gives a reliable consistency check on the quality of the 
data and the theory presented. 

It is evident, as shown in Fig. 2 and Table I, that espe- 
cially at T" > 30 K the data could not be accounted for by 
fewer than two incoherent oscillator modes superimposed 
on the electronic and coherent phonon contributions to 
the specific heat. The excellent correlation coefficient for 
the One-Debye-Two-Einstein mode scenario fit of the in- 
situ measurement {R = 0.99997) produced a very reliable 
fit, suggesting that in order to properly take account of 
that data we need a minimum of two "rattler" modes 
(see Table I). 
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FIG. 4. (color online). Close-up of the experimentally- 
measured density of states of Ybo.2Co4Sbi2 (the red line 
through the data points is a guide for the eye) exhibits a 
clearly-defined mode beyond the experimental resolution at 
a; « 5 meV, which is attributed to the dynamics of Yb, as 
shown in Fig. 1(b). 



The density of states G(uj) obtained from the inelas- 
tic neutron time-of- flight measurement on Ybo.2Co4Sbi2 
reveals the presence of low-lying phonons (w sa 5 — 25 
meV), separated from high-energy ones (w sa 27 — 35 
meV) (see Fig. 1). As shown in Fig. 1, there is an 
excellent agreement between features observed in our ex- 
perimental data (Fig. 1(a)) with those predicted by the 
density functional perturbation theoretical (DFPT) cal- 
culation on Ybo.25Co4Sbi2 (Fig. 1(b)). ^^ The modes ob- 
served between 27 meV and 35 meV are shown to be 
mainly due to vibrations of the Co atoms and arc in the 
same energy range as "iron modes" observed in the Fe- 
based skutterudite (Ce,La)Fe4Sbi2.^ 

A significant local maximum from the calculations of 
the phonon DOS is a low energy peak at w « 5.52 meV, 
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FIG. 5. (color online). Shown are close-ups of (a) 
the experimentally-measured G{u}) of Ybo.2Co4Sbi2 and (b) 
DFPT calculation of G(aj) of Ybo.25Co4Sbi2, both at higher 
transfer energies (13 - 19 meV), where a number of sharp fea- 
tures are resolved in the experimental data. The large blue 
and red arrows identify the locations of the peak [uj « 12.4 
meV) in the DOS of Ybo.2Co4Sbi2, which is identified with 
the second Einstein mode, Qe2- The small purple and pink 
arrows identify the locations of local modes at intermediate 
transfer energies {u ~ 10 - 22 meV) in the experimental data 
and DFPT calculation, respectively. 



which is attributed to vibrations of the Yb atoms (see 
Fig. 1(b)). A corresponding feature to the latter also 
appears in the measured DOS, as shown in Figs. 1(a) 
and 4, namely, a peak at 5.0 meV, which is resolved and 
well-separated from the band of Sb vibrations between 
7 meV and 23 meV (the high resolution of the CNCS 



instrument is readily visible in G{oj) vs. w close-ups in 
two separate tu ranges, as shown in Figs. 4 and 5(a)). 
The results from the fit with two Einstein oscillators of 
the in-situ specific heat measurement (see Figs. 2 and 3) 
corroborate the presence of a low-energy Einstein mode, 
Qei = 70 K [lljei = 6.0 meV), which corresponds to the 
measured and calculated values for the Yb peak at w w 
5 meV in G(aj) vs. uj (see Figs. 1 and 4). 

In previous heat-capacity studies on filled skutteru- 
dites a good consistency was observed between the weight 
of the first Einstein oscillator mode {A) with the theoret- 
ical prediction that each oscillator heavy element atom 
should contribute 1.47 J/K-mole of oscillator atoms. ^'^ 
In the latter papers the authors cite A = 1.21 J/K- 
mole in Lao.9Fe3CoSbi2.^'^ Ironically, in the case of 
Ybo.2Co4Sbi2, we would expect A to be 0.308 J/K-mole, 
whereas the fitted value is shown to be 1.62 J/K-mol - 
a difference by a factor of 5.25. The apparent similarity 
of A in Ybo.2Co4Sbi2 and Lao.9Fe3CoSbi2, despite the 
fact that the former has 4.5 times fewer rattler atoms 
than the latter, suggests that even the peak at 5 meV 
{Qei mode) could be a consequence of a local change 
in the Co4Sbi2 lattice. If one assumes that the nearest 
neighbors of Yb become loosely bound to the skutteru- 
dite lattice due to the filler, then their added degrees of 
freedom should account for the much greater A observed 
in Ybo.2Co4Sbi2. 

The pronounced appearance of the eigenmode at 5 
meV in Ybo.2Co4Sbi2 is attributed to the large neutron 
flux and improved resolution of the CNCS, but also to 
a large extent — to the large neutron scattering cross 
section of ytterbium (cr(Yb) = 23.4 bn, a/M = 0.13). 
The Yb mode is characterized by an energy range sim- 
ilar to other filler modes observed in a number of Ce, 
La and Tl-filled skutterudite systems*'^, and exhibits a 
striking resemblance to the results of a recent high reso- 
lution inelastic neutron scattering experiment and the- 
oretical DOS calculation on another Yb-filled system, 
YbFe4Sbi2.^'^ The latter system exhibits two separate 
filler modes, which manifest as distinct peaks at 4.9 meV 
and 5.7 meV in the experimental data, which are also 
seen in the theoretical calculations of Koza et al?^ . These 
are to be contrasted with our own results, which attest to 
a single peak in Ybo.2Co4Sbi2, both theoretically and ex- 
perimentally within resolution limits. We attribute this 
distinction between the guest modes of the two samples 
to differences in the curvatures of the Yb binding poten- 
tials in the FeSbs and CoSba cages, respectively. 

Furthermore, the specific heat measurement of 
Ybo.2Co4Sbi2 also seems to suggest the presence of an- 
other local mode at higher energy transfers (a; > 5 meV), 
beyond the Yb peak. In addition to 9bi, a second Ein- 
stein mode, Qe2 = 149.02 K, (w = 12.8 meV) is also 
accounted for from the functional fit through the specific 
heat data (see Fig. 2), and coincides with a peak in the 
DOS of Ybo.2Co4Sbi2 at 12.4 meV (see Fig. 5(a) and 
5(b)). The presence of Qe2 seems consistent with the 
two Einstein modes found in a number of other systems 



and has been related to the phonon glass. ^"'^ Upon closer 
inspection of the contributions of Oei a-nd Qe2 to the 
total Cp/T, 852 's part exceeds B^i's throughout the 
high temperature range, whereas below T « 27 K, O^i 
begins to dominate (see Fig. 3).^ However, the effective 
contributions of the two Einstein modes to the lattice 
thermal conductivity would be modified by the power- 
law frequency dependence of the attenuation of sound 
waves a {a"'""''' "'"'«=" oc w^ and a'°"ff**"*"°' '""'«=" o^uj^), 
as well as the absolute temperature, T, and the elastic 
constants of the system,^'' which could be obtained from 
a velocity of sound measurement on filled skutterudite 
systems. 

Another remarkable finding in our experiment is the 
presence of sharp modes {van Hove singularities) in the 
14 meV < w < 20 meV range (see Figs. 1(a) and 5(a)), 
which the DFPT calculation shows (Fig. 5(b)). These 
are mainly attributed to vibrations of Sb atoms (see Fig. 
1(b)). The 14 - 20 meV energy transfer range is distin- 
guished by 6 sharp peaks (shown in Fig. 5), and the one 
listed as " 1 " occurs at an energy transfer where Sb and 
Co peaks are coincident in the DOS of Ybo.2Co4Sbi2 (see 
Fig. 1(b) and Fig. 5). This suggests that peak "i" may 
be due to a hybridization of atomic degrees of freedom, 
as suggested by Vcithen and Ghosez and Feldman et al., 
whose calculations of the dispersion curves of CoSba ex- 
hibit a large number of nearly dispersionless optic modes 
that result in sharp features in the calculated phonon 
DOS."'^^'^^ What is surprising in our study is that we ob- 
served these sharp features, whereas other such experi- 
ments on CoSbs show only a broad peak in this energy 
range. ^ 

To further study the proposed guest-host coupled dy- 
namics of filled skutterudites, we also measured the 
neutron-weighted density of states of a related skutteru- 
dite prepared in the same way; Ce1.05Fe4Sb12.04- G{u) 
of Ce1.05Fe4Sb12.04 is shown in Fig. 6(a) and com- 
pares favorably with the previously-measured DOS of 
CeFe4Sbi2.^ A broad feature is present in the energy 
range 14 meV < w < 20 meV, and appears quite distinct 
from the sharp peaks observed in the Ybo.2Co4Sbi2 DOS 
(see Figs. 5(b) and 6(b)). Even though Cei.o5Fe4Sbi2 
seems to have a local maximum in the DOS centered 
around 14.3 meV (very faint bump), it appears to be 
less significant than the corresponding one (a peak at 
14.7 meV marked with a blue arrow in Fig. 5(a)) in the 
experimentally- measured DOS of Ybo.2Co4Sbi2. ujpeak is 
approximately 14 meV by the calculation and is marked 
with a red arrow in Fig. 5(b) (see Figs. 5(a), 5(b) and 
6(b)). On the other hand, the presence of ytterbium 
cannot be exclusively credited with the appearance of 
dispersionless modes in the intermediate energy range in 
the DOS of Ybo.2Co4Sbi2, since the density of states of 
YbFe4Sbi2, a different Yb-filled antimonide, shows a rel- 
atively broad feature in the same energy transfer range. ^'^ 

For the most part. Fig. 6(b) shows noisy data which 
indicate that either Ce1.05Fe4Sb12.04 does not exhibit dis- 
persionless modes or that even if they were to exist they 
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FIG. 6. (color online). (a) The density of states of 
Ce1.05Fe4Sb12.04 shows the same qualitative behavior as the 
results of a previous experimental and theoretical investiga- 
tion of CeFe4Sbi2.'-* (b) The DOS of Ce1.05Fe4Sb12.04 zoomed 
in on tlie same energy and amplitude scales as the correspond- 
ing DOS for Ybo.2C04Sbi2 failed to show well-defined disper- 
sionless peaks. The blue line is a guide for the eye. 



could not be resolved by the experiment. The broad 
noisy shoulder in CeFe4Sbi2.o4 bas also been observed 
previously by Koza et ai, with no additional features 
standing out.^ 

The DOS comprises the sum of the modes of the sys- 
tem over all reciprocal lattice vectors at a given energy 
exchange. Therefore, our experimental DOS represents 
a superposition of many different bands and the relative 



contributions of each of them. Thus, the important ques- 
tions are not only what is the scattering cross section of 
Yb versus Ce, since in the intermediate energy-transfer 
regime the theoretical calculations do not expect ytter- 
bium or cerium to contribute to the density-of-states, but 
rather ~ do dispersionless modes exist in the Ce-filled 
compound and how do they compare with the contribu- 
tions of coherent modes to the density-of-states. We at- 
tribute the noisy broad shoulder with no distinct peaks 
in Cei 05Fe4Sbi2.04 to a smaller relative amplitude of in- 
coherent modes compared to Ybo.2Co4Sbi2. 

The effects of the specific filler on the specific cage 
must be therefore considered. Ybo.2Co4Sbi2 seems to 
be a very special system, since it defies the traditional 
"rattling mode scenario" and seems to suggest a more 
complicated system than the ones previously described. 

The discrepancy between the theoretical calculation 
of spectral weight A and the experimentally measured 
one and the differences among densities-of-states between 
Ybo.2Co4Sbi2 and other similar systems suggest that a 
possible cause for the emergence of local modes in the 
10 - 20 meV energy transfer range are incoherent lattice 
vibrations of Yb nearest neighbors in the Ybo.2Co4Sbi2 
cage, since at those energies Yb is not expected to con- 
tribute directly. The DFPT calculation fully relaxes 
the filled structure, calculates the dynamic matrix at 
first and then performs the phonon spectrum calcula- 
tion. Thus, the DFPT calculations do include the effect 
of structure change, such as the structure relaxation and 
bonding strength change, due to the filling. In other 
words, these effects do contain the "coupling" of Yb rat- 
tling and local Sb vibrational modes in an implicit way. 
But no phonon-phonon coupling is explicitly included. 
In principle, there could be a change in the phonon DOS 
in the 14 - 20 meV regime as a result of Yb filling be- 
cause the Sb configuration of the cages does change in 
comparison with the unfilled systems. 

Altered lattice properties due to the presence of a 
heavy earth element could lead to localized "pocket" 
modes which manifest themselves as dispersionless modes 
at w = UJE2 = 12.4 meV, as well as in the 14 - 20 
meV regime, ^^'^^ i.e., low-energy modes appearing in the 
gap of a phonon-gapped system resulting from a pseudo- 
vacancy at an impurity site. Such modes result from the 
displacements of lighter atoms, such as Sb, which are 
nearest-neighbor atoms to the Yb filler. The vacancy 
which is left as a result of the displacement would be a 
mode associated with atoms moving in and out of these 
quasi-holes in the lattice. 

The appearances of a second Einstein mode (0_e2) and 
other sharp modes in the DOS of Ybo.2Co4Sbi2 in the 
14 meV < a; < 20 mcV range (where ytterbium is shown 
to not contribute to the phonon density-of-states) are at- 
tributed to altered properies of the skutterudite cage due 
to the addition of ytterbium. Given the high ZT of the 
filled skutterudites over their unfilled parents, the above- 
mentioned coincidence of dispersionless modes from the 
phonon density of states with a second Einstein mode in- 



dicates that the scattering of heat-carrying phonons oc- 
curs at a variety of energy scales and could also be owed 
to a novel mechanism which merits further research and 
investigation. 



IV. CONCLUSIONS AND SUGGESTIONS 

Our results have shown a pronounced peak in the 
phonon DOS of Ybo.2Co4Sbi2 at w ?» 5 meV, established 
from inelastic neutron time-of-flight scattering, specific 
heat measurements and theoretical calculations. This is 
interpreted as a " rattler" mode of Yb atoms vibrating in 
the cage of the CoSba host lattice. 

Novel contributions of this work include the obser- 
vations of the second Einstein oscillator in the DOS of 
Ybo.2Co4Sbi2 at uje2 = 12.4 meV, and the dispersionless 
modes in the w = 14 - 20 meV energy-transfer regime. 
Both are believed to arise as a consequence of ytterbium 
filling. It appears that incoherent modes from the cage 
itself are induced by the addition of Yb and these also 
participate in heat dissipation. 

Future acoustic experiments aimed at measuring differ- 
ences between the elastic constants of filled and parent 
cobalt antimonides as a function of frequency should be 
performed in order to understand the significance of each 
"rattler" mode to the total lattice thermal conductivity, 
kl-'^'^ In addition, future experiments aimed at measur- 
ing the DOS of filled cobalt antimonides as a function 
of filler concentration should try to assess the role of the 
corresponding filler on the number of intermediate en- 
ergy modes, their amplitudes and possible energy shifts 
in the latter. 

The further understanding of intermediate energy 
modes will most likely illuminate the importance of filler 
atoms in filled skutterudite thermoelectrics and allow the 
engineering of materials with higher dimensionless figures 
of merit. 
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